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Abstract 

In this paper we consider the question of sensor network 
coverage for a 2-dimensional domain. We seek to com- 
pute the probability that a set of sensors fails to cover 
given only non-metric, local (who is talking to whom) 
information and a probability distribution of failure of 
each node. This builds on the work of de Silva and 
Ghrist who analyzed this problem in the deterministic 
situation. We first show that a it is part of a slightly 
larger class of problems which is #P-complcte, and thus 
fast algorithms likely do not exist unless P=NP. We 
then give a deterministic algorithm which is feasible in 
the case of a small set of sensors, and give a dynamic al- 
gorithm for an arbitrary set of sensors failing over time 
which utilizes a new criterion for coverage based on the 
one proposed by de Silva and Ghrist. These algorithms 



and those that sought a non-deterministic approxima- 
tion but assumed significant capabilities of the sensors. 
For a survey of the literature, see 
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build on the theory of topological persistence 12 



The former approach requires a great deal of prior 
knowledge about the geometry of the domain and the 
exact location of the sensors. The latter, does not re- 
quire this exactness, but often requires a uniform dis- 
tribution of nodes or a high level of intelligence in the 
sensors. The main contribution of [9] was a criterion for 
coverage that requires none of these things. 

In the current paper, we take a middle ground and 
address the question of computing the probability of 
failure of the criterion of [9] given the probability of 
failure of each sensor. We show that a computing the 
probability of failure for a generalized set of complexes is 
NP-hard, but we give an algorithm which can be used to 
solve small instances of the problem, and an alternative, 
dynamic algorithm to give an early warning of potential 
failure. 
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1 Introduction 



The newly emerged field of Computational Topology 



12 continues to find ever increasing areas of applica- 
tion. Perhaps its most significant application so far has 
been in the use of topological data analysis (TDA) on 
a wide variety of datasets [12] [3] [6] , and has also been 
used effectively to find structure in image s [4] [TT] , shape 
in proteins and protein complexes [Tj [2] 15 and in 
many other areas. Recently, it was applied to sensor 
networks in 14 , 19] , 10, 16 and the current paper is 



an extension of [9] 

Topology enters the study of sensor network when we 
consider questions like coverage. When does a set of 
sensors effectively monitor a region and when are there 
gaps? Phrasing this geometrically, we start with a set 2 
of sensors x in a domain Act 2 where each can detect 
objects in a circular region of fixed radius r c , and we ask 
if the union of these discs covers all of A. This prob- 
lem has been studied quite a bit, but previous to [9j, 
most work fell into one of two groups - approaches that 
utilized geometric analysis to obtain an exact answer 
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Outline. 

The paper is organized as follows. Section [2] gives an in- 
troduction to Rips complexes and persistent homology, 
both of which will be used extensively in this paper. In 
Section 3 we summarize the problem and results of [9] . 
Section 4 adds the assumption that sensors have a prob- 
ability of failure, and Section [5] discusses the complexity 
issues of determining the probability that there is cov- 
erage of the domain. Section [6] presents a deterministic 
algorithm for those times when the set of sensors is small 
enough, and Section[7]gives a dynamic algorithm for use 
when the set of sensors is too large. 



Rips Complexes and Persistent 
Homology 



Let x be a set of points in K 2 and suppose that r > is 
given. We are interested in the topology of W r (x), the 
union of balls of radius r about the points of \- One 
can build a variety of complexes with vertices x that 
capture this topology, the simplest and most intuitive 
is the Cech complex which has a simplex < vo, . . . ,V}~ > 
whenever the balls of radius r/2 about the Vi have a 
non-trivial intersection. The nerve lemma tells us that 
Cech does indeed capture the topology of W r (x)> but it 
can be difficult to compute. In particular, it requires one 
to know the exact location of the points of x, a luxury 
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Figure 1: A visual representation of births and deaths 
of homology classes. The class a is born at K because 
it is not in the image of the map from It 
dies entering Kj because it was still not in the image 
of H{Ki_i) at H(Kj-i), but has merged with an older 
class upon entering H(Kj). 

that we do not have in this case. We must therefore use 
an approximation known as the Rips complex for the 
problem at hand. 

The Rips complex 7Z has a d-dimensional simplex 
cr whenever \\vi — Vj\\ < r for every pair of vertices 
Vi,Vj G a. Unfortunately, the Rips complex does not 
retain the homotopy type of W r (x), but what is lost in 
topological data is made up for in ease of computation. 
Furthermore, the only information necessary to build 
the Rips complex is the set of pairs of points whose 
distance is below a the prescribed threshold r. 

Since we will be considering nitrations of our simpli- 
cial complex and looking at how the topology changes 
with the change of the simplicial complex, a brief review 



of persistent homology is in order 12 



We begin with a filtration of a simplicial complex 1Z, 
given by a series of inclusions 

K 
II II 
Ko^Ki -^K 2 ^ *K n 

which induces maps on homology 

H{K) 



H(K ) — *-H{Kx) — *~H(K 2 



-H(K n ) 



We will use homology with 1*2 coefficients for the en- 
tirety of the paper. Consider how the sequence of 
homology groups changes as the simplicial complex 
changes. A class [a] G H p {Ki) is born at Ki if it is 
not in the image of the map H p (Ki_i) — > H p (Ki). 
This class dies entering Kj if once there it merges 
with an older class, i.e. lies in the image of the map 
HpiKi-t) H p (Kj). See Figure [I] 

To determine when classes are born and die in the 
filtration, we build the boundary matrix D. This is a 
square matrix with a row and column for each simplex 
in K, ordered with respect to the filtration (which en- 
sures that a simplex comes after all of its faces) . D is a 
0-1 matrix which has a 1 in location D[i, j] if and only 
if the ith simplex is a face of the jth. Applying the 
persistence algorithm as in |12| yields a reduced matrix 
R = DV , where D is the boundary matrix and V is an 



elementary matrix storing the column operations per- 
formed on D during the persistence algorithm. Here, a 
reduced matrix is one in which every column is either 
completely zero, or the lowest 1 in the column is not 
the same as the lowest 1 in any other column. We write 
Iowr(J) — i if the lowest 1 in column j of matrix R is 
in row i. 

To read the births and deaths from the matrix, note 
that a p-dimensional class is born with addition of a di- 
mension p simplex a if the column corresponding to a 
is completely zero. A representative for the class that is 
born is stored in the corresponding column of V. The 
class born at a dies with the addition of r if the lowest 
one of the column corresponding to r is in the row cor- 
responding to cr, in which case a and t are paired. If 
the addition of a simplex a gives birth to a class, it is 
called a positive simplex. Similarly, if the addition of a 
simplex r gives death to a class, it is called a negative 
simplex. 

The dc Silva-Ghrist criterion requires us to work with 
persistent relative homology. In this case:, we take a pair 
(7Z, J-), where T is a subcomplex of TZ, and a filtration 

T^TZ X ^TZ 2 -^ *1Z 

inducing maps on relative homology 
H{T, J=)-^H(Hi , T)-+H{n 2 ,F)-+ *H(1Z, T) 

and consider births and deaths in the usual way. This 
requires a slight modification of the boundary matrix by 
reordering the rows so that those simplices which are in 
the subspace J- are moved to the top of the matrix prior 
to performing the persistence algorithm. Then the ad- 
dition of a gives birth to a class if its column is cither 
zero, or the lowest one in its column corresponds to a 
simplex in T . Simplex r gives death to a class if its 
lowest one corresponds to a simplex which is not in J- ' . 
Other than these distinctions, computing persistent ho- 
mology in the relative case is the same as in the absolute 
case. 



3 The Coverage Criterion 

Working in a simply-connected domain in the plane, 
suppose that we have a set of sensors with a fixed ra- 
dius of coverage. Our goal is check that these sensors 
cover the whole domain. What makes this a challeng- 
ing problem is that we do not assume that we know the 
locations of the sensors. This means that standard ge- 
ometric techniques are not applicable. Instead we turn 
to topology to answer the coverage question by build- 
ing the Rips complex on the set of sensors, thought of 
as points in the plane. We can then use homology to 
check for holes in the coverage. 

Let x be the set of points corresponding to the loca- 
tion of the set of sensors in a compact connected domain 
A C K 2 which has a piecewise linear boundary. Sup- 
pose that each sensor has a fixed coverage radius r c > 0. 
The question is whether every point in A lies within 
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distance r c of some sensor in x- We do not use the 
distance r c to build the Rips complex, instead we add 
an additional capability to each sensor. Let r& > be 
fixed, with r& < \%r c for technical reasons. Each sensor 
is given a unique identification number to broadcast. If 
another node is within distance r^, it can hear the signal 
and identify the ID number, but it has no information 
about the location of the broadcaster. In particular, it 
does not know its direction or its exact distance, only 
that that distance is less than TV Whenever two sensors 
can hear each others' identification number, an edge is 
placed in the Rips complex. Higher dimensional sim- 
plices are then added when all of their faces are already 
there. 

The boundary of the domain A is taken to be piece- 
wise linear with a sensor at each of its vertices. The 
boundary is called the fence, and each node in the fence 
knows the identification number of its two fence neigh- 
bors, both of which are within distance 7v 

Summarizing, following [9], the assumptions are: 

1. Nodes x broadcast their unique ID numbers. Each 
node can detect the identity of any node within 
broadcast radius tv 

2. Nodes have radially symmetric covering domains of 
cover radius r c > r^j Vo. 

3. Nodes lie in a compact connected domain Acl 2 
whose boundary OA is connected and piecewise- 
linear with vertices marked fence nodes Xf- Non- 
fence nodes are called interior nodes, and denoted 

Xint' 

4. Fence nodes Xf are ordered cyclically and each v £ 
Xf knows the identities of its two neighbors on OA. 
These neighbors both lie within distance of v. 

Based on this information, build the Rips complex TZ 
with the fence J 7 as a subcomplex. With this setup, de 
Silva and Ghrist in [9] give their controlled boundary 
criterion for coverage: 

Theorem 1 (de Silva/Ghrist Criterion (dS-G)). // 
there is a nontrivial element of the relative homology 
group H2(TZ,J-) which maps to a nonzero class un- 
der the connecting homomorphism H2{TZ,J-) — > H\(!F), 
then the union of the disks of radius r c about the nodes 
contains all of A. 

The class [a] <G H2(TZ,J 7 ) is fundamental if it satis- 
fies the criterion of Theorem [l] but we stress that when 
there is such an element it is not necessarily unique. 
The term absolute cycle will be used for a class in 
H2(7Z,J-) that comes from Bi2(TZ), which is equivalent 
to saying that it maps to under the connecting homo- 
morphism. 

The assumption that r c > r^/VS is required to com- 
pensate for the fact that the Rips complex does not 
accurately reflect the topology of the cover. While this 
bound promises that holes in the cover appear also as 
holes in the Rips complex, we can still create examples 



where phantom holes appear in the Rips complex even 
though no hole exists in the cover itself. In a perfect 
world, this theory would be built on Cech complexes, 
however the lack of location data for the nodes makes 
this method impossible. 

4 Sensor Failure 

Over time, sensors have a likelihood of failure which 
increases the longer the system is in place, caused per- 
haps by malicious actions, environmental conditions or 
mechanical failure. As nodes fail, there are two possi- 
ble effects on the coverage: either the death of a subset 
of nodes creates a hole in the Rips complex, or the re- 
moval of the nodes does not affect the existence of a 
fundamental class. Once again, we emphasize that we 
are specifically not looking for the probability of failure 
of the cover over time, just the failure of the dS-G crite- 
rion. For this reason, we also assume that only interior 
nodes can fail. The loss of a fence node causes instant 
failure of the dS-G criterion, so there is nothing to check 
in this case. 

Let us start with a Rips complex pair (TZ, J-) built 
from a set of nodes x- At time t — 0, we assume we 
have a fundamental class [a] £ H 2 (TZ,J : ). If a set of 
interior sensors B C Xint fails, any simplex in TZ that 
has a vertex in the set B is lost. Therefore this subcom- 
plex, 7Zb , can be thought of as the largest subcomplex 
of 1Z that has x ~ B as its vertices. We could then de- 
termine whether TZ b fails the dS-G criterion by looking 
for a fundamental class in H 2 (TZb;J'), but this is a very 
narrow view of the problem. Much more information is 
available in a filtration that we will now construct. Note 
that it will contain TZb as one of its subcomplexes. 

Let \xint\ = n - Order the nodes so that B = 
{vi,- ■■ ,v k } and Xint - B = {v k+ i, ■ ■ ■ 7 v n }. Let Vi = 
{vi, ■ ■ ■ ,Vi} so that Vfc = B, and let TZi — TZvt be the 
maximal subcomplex of TZ with vertices X~Vi- Then 
the filtration 

T TZ 

ii ii 

TZ n — *~TZ n -\ — ^TZn-2 — **• • — *~TZ() 

induces maps on relative homology 

H(F, T)^H(TZ n - U f)^H{TZ n -2,T)— -H{H, F) 

An example of this filtration is illustrated in Figure [2] 
Intuitively, we expect that discovering a fundamental 
element at any point in this sequence implies that there 
is a fundamental element in any subsequent group. 

Lemma 1. Let A c B be subsets of Xint- Then if 
[f3] G H2 (TZb, J 7 ) is fundamental, its image under the 
map 

H 2 (n B ,F) H 2 (n A ,F) 

is also fundamental. 
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Figure 2: A small scale example of the Rips complex TZ 
built from a set of points in the plane. The outer ring 
of nodes labeled with letters is the fence J ' . 



Proof. Consider the commutative diagram 



H 2 (n Bl T) 

8 



9 



(1) 



where the horizontal maps are induced by the inclu- 
sion TZb C TZ a and the vertical maps are the bound- 
ary maps. Since [/3] 6 H 2 (JZbt^F) is fundamental, 
d[/3] 7^ in Hi (J 7 ). Since the diagram commutes, 
9i*[/3] = d[/3] and hence is nonzero. This also implies 
that i*[/3] is nonzero, so it is a fundamental element of 

H 2 {n Al F). □ 

This lemma shows that if TZb passes the dS-G crite- 
rion, then TZa passes the dS-G criterion for all A C B. 
But it also shows that if TZb fails the dS-G criterion, 
then TZa fails the dS-G criterion for all Ad B. Thus if 
B is a death set, and B C A, Lemma [l] implies that A 
is also a death set. This leads us to make the following 
definitions: 

Definition 1. A set B C Xint whose removal causes 
failure of the dS-G criterion is called a death set. A 
death set B C Xint * s a minimal death set if no subset 
of B is itself a death set. 



Definition 2. // the removal of B does not cause fail- 
ure, we call B a cake set. A cake set B C Xint is a 
maximal cake set if no superset of B is also a cake 
set. 



As we will show in section |6.2| minimal death sets 
are directly related to the failure of the dS-G criterion. 
However, we first look at the issues arising from the 
complexity of the problem. 

5 Complexity Issues 

The first issue to address is whether this problem is com- 
putationally complex. In this section, we will show that 
in fact it is part of a slightly larger group of problems 
which are NP-hard, more specifically #P-complcte. 

5.1 Use of the 2-skeleton 

A simplifying step is to work with the 2-skeleton TZ 2 of TZ 
rather than the full complex. To justify this, notice that 
passing to the 2-skeleton does not affect our observance 
of the dS-G criterion: 

Lemma 2. The dS-G criterion is satisfied for TZ if and 
only if it is satisfied for TZ 2 , where TZ 2 is the 2-skeleton. 

Proof. Consider the following diagram built from the 
long exact sequences for the pairs (TZ,T) and (TZ 2 ,^) 
and the maps induced by the inclusion TZ 2 — > TZ: 



H 2 {TZ 2 ,F) 



H 2 {TZ,F) 



H^TZ 2 ) 



Y 

H^TZ) 



If the dS-G criterion is satisfied for TZ 2 , then there is an 
a G H 2 [TZ 2 , F) such that d(a) is nonzero. Clearly i*(a) 
also satisfies the dS-G criterion since 9(i*(a)) = da. 

Now assume that the dS-G criterion is satisfied for TZ. 
Then there is a /3 £ H^TZ,?) such that df3 is nonzero 
in Hi(T). Here it is important to note that since TZ 2 is 
the 2-skeleton of TZ, H^TZ) = H 1 (TZ 2 ). As the top and 
bottom rows are exact with the last two groups equal, 
and since df3 £ H\{F) maps to zero in Hi(TZ), it follows 
that it also maps to zero in Hi(TZ 2 ). Because the top 
row is exact, there is an a € H 2 (TZ 2 ,J-) which maps to 
dp, and hence satisfies the dS-G criterion. □ 

This lemma implies that the sets of death sets, min- 
imal death sets, cake sets, and maximal cake sets are 
equivalent to their counterparts when computed in the 
2-skeleton. It also implies that the probability of fail- 
ure of the dS-G criterion in the full Rips complex and 
the probability of failure of the dS-G criterion in the 2- 
skeleton are the same. And lastly, in TZ 2 there is exactly 
one cycle representing each homology class (Z 2 = H 2 ). 
For these reasons we will simplify notation and write TZ 
for the 2-skeleton of the Rips complex for the remainder 
of the paper. 
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5.2 #P-Complete 

The class of problems defined as #P-complete was in- 
troduced by Valiant in [17]; they form a specific class of 
NP-hard problems. Typically, #P-complete problems 
are concerned with counting how many of something 
exists whereas general NP problems just ask if some- 
thing exists. Problems which are #P-complete likely do 
not have polynomial time algorithms. 

To show that a problem is #P-complete, we reduce 
one difficult problem to another. Reducing problem A 
to problem B means that we take any instance of prob- 
lem A, use it to create an instance of problem B, and 
conclude that the answer to solving problem B gives 
an answer to problem A. To prove NP-completeness 
or #P-completeness, both turning an instance of prob- 
lem A into an instance of problem B and returning the 
answer to problem A given the solution to problem B 
must be done in polynomial time. If we can reduce A 
to B in polynomial time, we write A <p B. 

A reduction from A to B is called parsimonious if the 
number of solutions for A is in one-to-one correspon- 
dence to solutions for B. This is an important property 
for proving that problems are #P-complete since we 
need to be able to count the number of solutions of A 
based on the number of solutions of B. 

In order to show that our sensor network problem is 
$F-complete, we need to find a polynomial time, parsi- 
monious reduction from a #P-complete problem to our 
problem. 

It has been known for several decades that the com- 
puter science problem of network reliability is #P- 
complete j8 13 . We will specifically work with the two 
terminal network reliability problem as defined in |13| . 
An instance of the problem is a graph G — (V, E) with 
marked vertices {s,t}, a rational failure probability p e , 
< p e < 1 for each edge e € E, and a positive rational 
number q < 1. Then, assuming edge failures are inde- 
pendent of one another, we ask whether the probability 
that s and t have a path with no failed edge is greater 
than or equal to q. The fact that this problem is hard 
in the class of counting problems comes from needing to 
count the possible paths from s to t when determining 
the probability of failure. 

This problem has striking similarities to ours, and the 
closeness is even more pronounced when we look at it 
in the following way. Considering the graph G as a one- 
dimensional simplicial complex, a path in G with end- 
points at s and t is a fundamental class in H\(G 1 {s, £}), 
where fundamental means that the boundary of the 
class is homologous to [s] + [t] in Hq({s, t}). 

Our goal is to reduce network reliability to our prob- 
lem, which we will therefore call 2-dimensional network 
reliability. An instance of the problem is a simplicial 
complex X with a subcomplex Y that is homeomor- 
phic to S 1 . We also have rational probabilities of fail- 
ure p v , < p v < 1, on the vertices not in Y, and a 
value < q < 1. We ask the following question: Given 
the fact that failures of vertices are independent of each 
other, is the probability that we have a fundamental class 



a e H 2 (X,Y) at least q? 

Notice that our definition of the problem takes no ac- 
count of the geometry inherent in the originally defined 
problem as we are ignoring the fact that we obtained 
this simplicial complex from a set of points in M 2 , and 
thus we are proving a larger class of problems to be 
#P-complete. To prove that 2-dimensional network re- 
liability is #P-complete, we must take an instance of 
the 1-dimensional network reliability problem, turn it 
into an instance of the 2-dimensional case in polyno- 
mial time, take the solution given there and turn it into 
an answer to the 1-dimcnsional case in polynomial time. 

Theorem 2. 2-dimensional network reliability is #P- 
complete. 

Proof. Consider a finite graph G with vertex set V and 
probability of failure p e given on each edge e. We will 
construct a 2-dimensional simplicial complex X with a 
subcomplex Y = S 1 so that there is a one-to-one corre- 
spondence between paths from s to t in G and funda- 
mental classes of H2(X,Y). This correspondence will 
also preserve the probability of failure of the class, so 
this will imply that the probability of failure in the 1- 
dimensional case can be computed by determining the 
probability of failure in the 2-dimensional case. 

Suppose that G has n vertices. Order these vertices 
by choosing a map r : V — > K that sends each vertex 
to a distinct integer in {1, • • • , n}, with r(s) = 1 and 
r(t) = n. Extend r to all of G by linear interpolation 
over each edge, and subdivide G by adding vertices at 
all points of r -1 ({l, •■• ,n}). Call the result G", the 
map r' : G' — > is now piecewise linear. If an edge e 
of G is subdivided into k subedges in G', we set the 
probability of failure for one of the subedges equal to p e 
and the rest equal to 0. See Figure [3] for an example 
of building this graph. 

Form the complex 

G'x//~ 

where (x,z) <~ {x',z') iff r(x) = r(x') and either 
z = z' = 1 or z = z' = 0. Note that this collapses 
the top and bottom graphs each onto a separate copy of 
the interval [r(s),r(i)]. To make this a true simplicial 
complex, divide each rectangle of the form ex I into tri- 
angles by placing a vertex at the barycenter and adding 
the obvious four new edges and four new triangles. The 
resulting complex will be called X. Then define Y to 
be the subcomplex G' x {0, 1} / ~ together with the 
two edges s x I and t x I; Y is homeomorphic to S 1 by 
construction. (In Figure [3] we have not subdivided the 
rectangles to keep the picture uncluttered.) 

Set the probability of failure of vertices that were 
added to the centers of the rectangles equal to the prob- 
ability of failure of the edge of G' from which they arose. 
Notice that failure of one of these vertices leads to re- 
moval of the interior of the corresponding rectangle. 

It is obvious that each path in G gives rise to a funda- 
mental class in H2(X,Y). To prove the opposite, first 
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r(u) = 2 




r(v) = 3 



r(u) = 2 r(b) = 3 




r{a) = 2 r(v) = 3 




Figure 3: G is the instance of the 1-dimensional Net- 
work Reliability problem with an ordering r placed on 
the vertices. This map can be extended by linear inter- 
polation over each edge. We then subdivide the edges 
at all points of r _1 ({l, • • • , n}). Finally, we define the 
complex X = G X 1 1 ~ where (x,z) <~ (x',z r ) iff 
r(x) = r(x') and either z — z' = 1 or z = z' = 0. 

recall that since X has no 3-simplices, each fundamen- 
tal class in H 2 (X, Y) has a unique representative cy- 
cle (B 2 (X,Y) = 0). Furthermore, since we are using 
homology with Z 2 coefficients, each class is simply a 
subset of the set of 2-simplices in X. A relative cy- 
cle a £ Z 2 (X,Y) has the added property that an even 
number of 2 simplices in a contain any edge of X — Y 
and a fundamental class must have da equal to the sum 
of all of the simplices of Y. Notice also that if a con- 
tains any 2-simplex from a rectangle, it must contain 
all four 2-simplices from that rectangle, so it is equiva- 
lent to think of a as a set of rectangles from before the 
subdivision. 

Our conclusion now follows easily. Every rectangle in 
a determines a unique edge of G' . Since T is covered 
exactly once, no two edges of G' that are equivalent 
under ~ can occur as edges of rectangles in a. Since a is 
a relative cycle, each vertical edge must lie on either or 
2 rectangles, so the edges patch together to give a path 
from s to t. Hence there is a one to one correspondence 
between the two sets. 

Since we set up each rectangle to have an equal prob- 
ability to that of its corresponding edge, the probability 
that a fundamental class is still in H 2 {X, Y) is equal to 



the probability that the corresponding path in G is still 
functioning. Thus, if we could compute the probability 
of failure in X in a reasonable time frame, the solution 
would give the probability of failure in G in a reasonable 
time frame. Since the latter problem is #P-complete, 
our 2-dimensional version is also #P-complete. 

□ 

6 A Deterministic Algorithm 

Now that we know that the general problem is #P- 
complete, we strive to find ways to work around the 
computational complexity issues. We will first show 
that, given a set of sensors which is relatively small or 
at the very least relatively sparse in the domain, we can 
write a deterministic algorithm to compute the prob- 
ability of failure of the system. In section [7J we will 
consider the modified problem of predicting failure as 
sensors in the system fail. 

6.1 The Hasse Diagram 

Consider a set of sensors x hi A. Recall that edges are 
added when sensors are within of each other, and 
2-simplices are added wherever all three edges have al- 
ready been included. 

Consider all possible subsets A C Xint an d as before 
construct the Rips complex TZa, the largest subcom- 
plex of 1Z which does not utilize the nodes in A. The 
collection of these Rips complexes forms a poset un- 
der inclusion, where A C B gives the reverse inclusion 
TZb C TZa- Arrange all of these Rips complexes into a 
Hasse diagram, as shown in Figure [4] for the example in 
Figure [2j Here we place TZa in the row indexed by the 
number of elements in A, and we have shaded all the 
complexes TZa which fail. A line is drawn between TZa 
and TZb if A is obtained from B by removing a vertex. 



n 




Tli K 2 Tl-i Ki 
TZ-12 TZ\i TZia 7^23 TZ 2 a ^34 




TZ\2iA 



Figure 4: Hasse Diagram. A Rips complex TZa is placed 
on a row according to the size of A and lines are drawn 
to show inclusion between complex is neighboring rows. 

In Figure [4] notice that, if TZb fails the dS-G criterion, 
the Rips complexes for all its supersets of B do as well, 
so all its successors are also shaded. This means that 
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when searching for failures we do not have to check every 
possible subset. Using breadth first search from 1Z, we 
only need check complexes where all the predecessors 
are cake sets. From this search pattern, if it is necessary 
to check the set in the first place then it is not just a 
death set but a minimal death set. This means that 
there is no post processing needed to determine the list 
of minimal death sets. 

Given this setup, we now consider the probability of 
failure of the dS-G criterion. 

6.2 Probability of Failure 

Let Xi be a random variable which gives the time of 
death of node Vi. In many cases, Xi will be an expo- 
nential random variable, but this has no effect on our 
result so we make no such assumption. We do however 
assume that Xi and Xj are independent for i ^ j. 

Let Sa be the random variable which gives the first 
time at which all nodes in the set have failed, clearly 
Sa = max{Xi,X2, • • ■ ,Xk}. Because the failures of 
the nodes are independent events, we have 

P(S A < t) = F(max{X u X 2 , ••• ,X k ] < t) 

= ¥{X 1 < t)P(X 2 <*)••• P(X k < t). 

Next, let V = {Ax, • • • , A r } be the collection of all death 
sets, not necessarily minimal. Let C be the random 
variable which gives the time of failure of the dS-G cri- 
terion for the system. The value of C gives the first 
time that all of the nodes in one of the Ai have failed, 
i.e. C = min{SA 1 , • • • ,Sa t }- Hence the probability 
that the system has failed the dS-G criterion by time t 
is given by 

P(C<t)=P(mm{S Al ,--- ,S Ar }<t). 

Unfortunately, these events are not independent since 
many of the death sets have non-trivial intersections. 
On the other hand, we can organize T> using the concept 
of death chains. 

A death chain is a sequence Ax C A2 C ■ ■ ■ C A q of 
death sets Ai € T>. A maximal death chain is a death 
chain which cannot be increased in length, either by 
inserting any intermediate set between elements of the 
chain, or adding any sets to either end of the chain. Note 
that such a chain starts with a minimal death set, ends 
with A n — Xint j and the number of elements increases 
by exactly one going from Ai to Ai+x . If we read off the 
Hasse diagram as in Figure |4j we see that a maximal 
chain is a path which goes from a minimal death set to 
the bottom of the diagram. 

We will call a maximal death chain with A as the 
minimal element an A-chain. This allows considera- 
tion of all possible maximal chains from the complex 1Z, 
grouped by first element: in our case yl^-chains through 
A^-chains. Using this organization of the death sets 
gives the following theorem: 

Theorem 3. Let {Ax,-- - ,^4d} be the set of minimal 
death sets for the Rips complex {TZ,T). Then the prob- 
ability that the complex has failed by time t is equal to 



P(Failure by time t) = P (min{,SUJ < tj . 

Proof. If T> is the set of all death sets, some set B in D 
must have failed in order to cause failure of the dS-G 
criterion. Every set in T> contains one of the minimal 
death sets Ai, so every set in T> is in at least one Ai- 
chain. 

It is obvious that given any Ai-ch&in Ch, where the 
time at which any of the death sets in the chain have 
failed is given by Y, satisfies 

F(Y<t) = ¥{S A < t). 

From this we can also conclude that any set of A-chains 
{Ch 1 , • • • , Ch fe }, where each respective time of failure is 
given by Y 1 , ■ ■ ■ , Y k , satisfies 

P =P(SA<t). 

Therefore, instead of asking for failure of the dS-G 
criterion, we can then ask for the time when at least 
one of the sets in at least one of the chains has failed. 
Hence 

P(Failure by time t) = P ^{5% < i}J 

= P {l)i S ^ ^ *}J 
= p(min{SUJ < tj . 

□ 

This means that the probability of failure of the sys- 
tem can be computed given the minimal death sets. No- 
tice that there is still work to be done since the mini- 
mal death sets may have intersections. However, since 
we assume that we have a small number of sensors that 
are well distributed and are not extremely dense in the 
domain, the number of intersections will be small and 
therefore this later computation is feasible. Thus, we 
seek an algorithm to compute the minimal death sets 
although from our knowledge that the problem is #P- 
complete, we expect that this algorithm will be expo- 
nential in the worst case. 

6.3 Death Sets Algorithm 

When constructing an algorithm to determine the set 
of minimal death sets, the search space is the Hasse 
diagram described in Section |6.1| Since this has size 
2lx<»t| ) we expect this is the source of our complexity 
issues. 

Search the Hasse diagram using breadth first search. 
This will exploit the property that if B is a death set 
and B C A then A is also a death set. Hence, to find 
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the minimal death sets, we must only check complexes 
where every predecessor is still a cake set. If we are 
forced to check all of the nodes of the Hasse diagram, 
then we will need to check 2l Xi71t l complexes. However, 
with our assumption that we do not have a very dense 
set of sensors, it should take removal of a small set of 
sensors in order to break the dS-G criterion. This means 
the size of minimal death sets will be relatively small, 
and thus they will be close to the top of the Hasse di- 
agram. More importantly, it means there will be rel- 
atively few of them so out output size will not be too 
large. 

Given this method to work through the sets, we need 
an efficient way to check the dS-G criterion. Consider 
the subcomplex TZa , thought of as the point in a filtra- 
tion of TZ where all simplices have been added except 
those which have vertices in A. Order the simplices so 
that those in the fence come first in the ordering. Ini- 
tial intuition says that in order to talk about failure of 
the dS-G criterion when the nodes in A are removed, 
we should filter TZ so that all nodes, edges, and trian- 
gles which have any vertex in A are last in the ordering. 
We could then construct the boundary matrix for this 
ordering, cut off the final columns corresponding to sim- 
plices which would be gone if the vertices in A failed, 
and reduce the resulting matrix in order to read off the 
homology of (TZa,J~). 

However, this turns out to be much more work than 
is needed. If we add all degree 1 and degree simplices 
from the start, even if they have interior nodes which 
we assume to have failed, the failure of the dS-G crite- 
rion is not affected. The following expresses this and is 
elementary to prove. 

Lemma 3. Let Y be a 1- dimensional simplicial com- 
plex, X a 2- dimensional simplicial complex whose vertex 
set may intersect nontrivially with Y , and J- C X , then 
H 2 (X,F)^H 2 (XUY,F). 

This implies that the time in the filtration when we 
add the 1-simplices is irrelevant to the homology group 
we are interested in, namely H 2 (TZ,J-). Thus, we can 
order our filtration so that all the 2-simplices are at the 
end and consider the failure of a node as the failure only 
of the 2-simplices which contain it as a face. 

Given this filtration, we reduce the matrix D via the 
persistence algorithm (see, e.g., [12]), and consider the 
rightmost columns, which correspond to 2-simplices. If 
the row for a simplex has no lowest one in a row below 
those corresponding to the fence simplices, the addition 
of that 2-simplex creates a new class in H 2 {TZ,F). If 
this column is not entirely and has a lowest 1 in a 
row corresponding to a fence simplex, that class has a 
boundary which is nonzero in H\(J-). Thus, our dS-G 
criterion reduces to looking for a column corresponding 
to a 2-simplex which has a lowest one in a row corre- 
sponding to a fence simplex. 

This shows that we can quickly determine whether a 
complex satisfies the dS-G criterion once we have de- 
termined the correct filtration for TZa and have reduced 



the matrix D. We would like to not have to rewrite 
and re-reduce the matrix R for each complex TZa to be 
checked. So, let us determine an efficient way to swap 
all the 2-simplices that have a vertex in our failure set A 
to the end. For this, we turn to [7], which gives an algo- 
rithm to quickly update and maintain the properties of 
R and U, where D = RU, as we swap columns (Notice 
that U — V^ 1 from the earlier discussion in Section |2j. 
This means that we do not have to rerun the persistence 
algorithm each time to reduce the matrix D. 

Let D = RU be an i??7-decomposition. That is, 
R is the reduced matrix, and U is upper triangu- 
lar. Let P be the matrix which swaps rows i and 
i + 1, so that PDF is the boundary matrix with sim- 
plices Ui and (Tj+i switched. This can be written as 
PDF = (PRP) (PUP), so we need to determine when 
PRP is not reduced and PUP is not upper triangular. 
Notice that PRP is not reduced if and only if there are 
columns k and I with low^(fc) = i, Iowr(1) = i + 1, and 
= 1. PUP is not upper triangular if and only if 
U[i,i + 1} = 1. 

In [7j , cases are split into whether Oi and <jj are posi- 
tive or negative. The one case that will not occur for us 
is the possibility that Ci and ai+ \ are positive, and there 
are rows k and I with low#(fc) = i and \owr(1) = i + 1. 
This would imply that 01 and <Jk are 3-simplices whose 
additions kill the classes born by the addition of cr, and 
<Ji+i- As we are assuming TZ is a 2-dimensional sim- 
plicial complex, this case is impossible, hence we can 
disregard it. With respect to the other cases, we can ei- 
ther swap rows and columns in R and U, hence replace 
them with PRP and PUP with no issues, or we must 
replace them with PRWPW and WPWUP, where W 
is the matrix which adds column i to column i + 1 in R. 
We replace R and U with PRP and PUP if 

<C> Ui and <7i + i are both positive simplices, 

<C> Ui and Ui+i are both negative simplices and 

U[i,i+1] = 0, 

<C> <jj is negative, cr-j+i is positive, and 
U[i,i + 1] = 0, 

()> Oi is positive and <7i+i is negative. 

We instead replace R with PRWPW and U with 
WPWUP if 

<0> Ui and are both negative and 
U[i,i + 1] = 1 

<) cr,; is negative, cr^+i is positive, and 
U[i,i+1] = 1. 

Hence, the columns of 2-simplices which correspond to 
failure of specific vertices can be quickly swapped to the 
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Fence Simplices Remaining 1-skeleton 




Column re 



of fundamental class 



resenting birth 



R 



Figure 5: The column representing the birth of a fun- 
damental class. It corresponds to a 2-simplex and its 
lowest 1 corresponds to a simplex in J- . 



end of our matrix. Now that matrix R corresponding 
to the filtration placing the simplices in A at the end 
has been reduced, we can check whether our complex 
TZa passes the dS-G criterion. In the language of per- 
sistence homology, the dS-G criterion is looking for a 
column which represents a positive simplex, and there- 
fore a simplex which adds a new class to H2{7Za, J 7 ). 
Additionally, the boundary of this new class is nonzero 
in Hi(!F). In terms of the matrices, we need to find a 
column i in the reduced matrix R which is has a lowest 
1 corresponding to a simplex in J- . A representation of 
this column is in Figure [5| Our algorithm is as follows. 



CakeOrDeath(ft, T) 

Given: The boundary matrix D with filtration 

{Fence, Remaining 1-simplices, 2-simplices } 
Reduce RU = D 

for A EXint, in the order of BFS in the Hasse diagram: 
Swap all columns corresponding to 2-simplices with 

a vertex in A to end of matrix D, and maintain 

R' and U'. 
if there is not a column as in Figure [5j 

mark A as 'Minimal Death.' 
else: 

mark A as 'Cake.' 
endif 
endfor 

6.4 Complexity of Algorithm 

As a beginning aside, we point out why we chose breadth 
instead of depth first search. BFS has the property that 
we will only ever check sets which are cake or minimal 



death. On the other hand, DFS would require post- 
processing to determine which of the death sets found 
were minimal death sets. The perk of DFS, however, 
is that it requires less matrix swaps since multiple sets 
A can be labeled as cake or death by reading off of one 
matrix. As we wish to have less post-processing, we 
choose to use BFS for our algorithm. 

Assume that the matrix D is stored as a 
sparse matrix. This is done with a linear array 
of lists D[l,"- , to] where m is the total num- 
ber of simplices in the 2-dimensional complex 
1Z. Also, assume the ordering of the simplices is 
{fence, vertices and remaining 1-simplices, 2-simplices}. 
Each entry D[i] in this array stores a linked list 
denoting the locations of the codimcnsion-1 faces of Cj, 
or equivalently, the Is in column i of the full matrix D. 
This not only speeds up the operations required on the 
matrix, but reduces the storage size of D to 0{m). 

There are four major parts to the algorithm. The first 
is to reduce D = RU. Using the persistence algorithm, 
this takes time at most 0(m 3 ) 



See 
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for details. 



For each complex to be checked, all necessary columns 
must be swapped to the right side of matrix R. If there 
are t 2-simplices, at worst there are t 2 /A swaps to be per- 
formed. While each swap has at worst an 0(m) running 
time, from |7] there is an amortized time proportional 
to the number of Is in the affected rows and columns, 
so this is 0(1). This step therefore has an amortized 
cost of 0{t 2 /A). 

Next, we check for a fundamental class. If this is 
done in TV complexes, then N < 2 Xint ~ 1 . (Recall that 
because we assume that the dense set of sensors is not 
dense, N will likely be much smaller than 2l Xi ™*l.) 

If a vector Low giving the location of the lowest 1 
in each column is maintained throughout the process 
of swapping, easily done via the cases in [7], it only 
takes 0(t) time to check for a column which fits our 
requirements. Hence for each death set, the amortized 
running time is 0{t 2 ), and so the overall running time 
is 0(m 3 + Nt 2 ). 

To make this running time feasible, one needs to keep 
N under control. The easiest way to do this is to have 
a sparse set of sensors. For example, suppose the area 
of the domain A is R and we have n sensors. Each 
sensor covers an area of irr 2 , so the sensors cover a total 
area (double counting overlap) of nirr 2 . This means 
the expected number of sensors covering any point is 
x = nirr 2 /R. Therefore the death sets should be of 
size approximately x, so it is only necessary to check 
complexes through about the x th row. Thus N ss (™) + 

In conclusion, computing the probability exactly will 
be easier if the set of sensors is sparse, but what is gained 
in exactness of the computation is lost in the robustness 
of the network. 
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7 A Dynamic Algorithm for a 
Monitored System 



Consider the diagram 



h 2 (tz w ,t) ^ H X (F) ^ H X (TZ W ) 



Suppose we are in a situation where the deterministic 
algorithm is not feasible. Computing the probability 
of failure exactly is an NP-hard problem, and thus an 
exact computation is essentially impossible when the 
set of sensors is large. Instead, assume that a central 
monitoring station receives information as to whether or 
not each sensor has failed. In this case, a more practical 
question is to ask when the system is getting close to 
failure and so we seek a dynamic algorithm to predict 
which nodes would cause failure of the criterion should 
they fail soon. To do this, we will create a new criterion 
built from the old which will give an early warning for 
failure. It essentially gives a flag on each interior vertex 
warning that its failure would probably cause failure of 
the dS-G criterion. 

For technical reasons, we will assume in this section 
that the domain is convex. A domain that is not convex 
can have a radius where the Rips complex has a non- 
trivial class in H X (TZ) even though (TZ, J 7 ) passes the 
dS-G criterion. This assumption is much stronger than 
is necessary since all we really need is that Hi(7Z) = 
whenever (TZ, JF) passes the dS-G criterion. 

The main idea for the new criterion comes from the 
following theorem. 

Theorem 4. Assume that the pair of simplicial com- 
plexes (TZ,F) passes the dS-G criterion, H 2 (TZ) = 0, 
and w is a vertex in TZ — T . Then (TZ W ,F) passes the 
dS-G criterion if and only if Hx(Lk(w)) = 0. 

Proof. Assume that (7Z W , JF) passes the dS-G criterion 
and that the domain is convex, which implies that 
H\{TZ W ) is trivial. Mayer- Vietoris for TZ — 1Z W U St(iy) 
gives the exact sequence 

H 2 (1Z) ^ ITi(Lk(ti;)) ^ H X {TZ W ) 8 Hi(St(w,ft)) 

Thus, since H 2 (TZ), H X (1Z W ) and Hx(St(w,7Z)) are all 
trivial, #i(Lk(u>)) = 0. 

Assume on the other hand that Hi(Lk(w)) = 0. Note 
that if (TZ,F) passes the dS-G criterion, Ho(Lk(w)) = 0. 
Using the long exact sequence of the pair (TZ W , Lk(w)), 
we have 

>Hx{1Z w ) >Hx{1Z w ,Lk(w,1Z)) ^0 

hence the middle map is an isomorphism. 

By excision, Hi(1Z w ,Ul(w,1Z)) = Hx(K, St(w,H)). 
Using the knowledge that St(w,7Z) has the homotopy 
type of a point for the first isomorphism and the long 
exact sequence of the pair (TZ, •) for the second, we have 

Hi(K, Si(w,1Z)) = Hx(1Z, ■) = Hi(K). 

Thus, Hx{TZ w ) Hx(TZ). 



I V I 
H 2 (TZ, F) Hx (F) > H t (TZ) 

Since (TZ, IF) passes the dS-G criterion, there is a fun- 
damental class a in H 2 (TZ, J 7 ). By definition, it maps to 
da which is nonzero in Hx(F), and since the two rows 
are exact, da maps to in H X (TZ). The last two verti- 
cal maps are isomorphisms, so there must be a nonzero 
(3 in H 2 (TZ W , J 7 ) which maps to da under the top left 
horizontal map, and therefore (TZ W ,J-) passes the dS-G 
criterion. 

□ 

Notice that the assumption that H 2 (TZ) — is only 
used for one direction of the theorem: if (TZ W , J-) passes, 
then we have a link with trivial first homology. De- 
spite being counterintuitive, it is possible for a set of 
points in the plane to have a non-trivial second homol- 
ogy group [5] . The expectation is that this event will not 
be frequent, but it must be kept in mind as we create a 
new criterion in this monitored set up. 

7.1 The new criterion and complexity 

Given Theorem [4j we propose a new criterion to com- 
plement the de Silva - Christ criterion: 

Definition 3 (Link Condition). If an interior vertex w 
has H\(Lk(w)) ^ 0, we say it is flagged. Otherwise, we 
say it is not flagged. 

The idea is that if a vertex is flagged, there is a chance 
its removal will cause failure of the dS-G criterion. If it 
is not flagged, then its removal can do no harm. 

With this definition in mind, we give a dynamic al- 
gorithm to follow as nodes fail: 

MonitoredSystemFailure(7?., J-) 

Given: Simplicial complex pair (TZ, J-) 

Check dS-G criterion (We assume that this initial 

check will always pass) 
Compute link of each vertex w and Hx(Lk(w)) 
if Hx(Lk(w)) = 0: 

Mark w as flagged, 
endif 

if vertex v fails: 
if v is flagged: 

Update matrix R to remove dead simplices 

Check dS-G criterion 

if (TZ W ,F) fails the dS-G criterion: 

Break 
endif 
endif 
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Update links of vertices 
Compute Hi(Lk(w)) for w whose link has 
changed 

Mark or unmark w as flagged according to Hi 
computation, 
endif 



This algorithm turns out to be polynomial in the 
number of simplices in the worst case. We will split 
the complexity computation into two parts: the initial- 
ization step, done before any vertex has failed, and the 
time taken for the algorithm for each failed vertex. 

Initialization 

Let m be the number of two simplices of dimension < 
2. In section 5.1 we showed that this is the highest 



< a, v >G X w . Therefore, a € ~Lk(v,X w ), and equation 
12 follows. □ 

This lemma implies that the only update needed af- 
ter the failure of a vertex is to delete any simplices in 
the link which were also deleted in the simplicial com- 
plex. In the worst case, the link of every vertex must be 
updated, the first homology recomputed, and the flag 
remarked as needed. Since the size of each link is at 
most m, this step takes 0(n(m + m 3 )) = 0(nm 3 ). 

If this sequence of events happens for every n, this 
second part of the algorithm takes time 

O ( n(t 2 + t + mn + m 3 n) ) = 0(m 3 n 2 ). 



dimension simplex needed for the dS-G criterion, and 
since the link condition only looks at the first dimension 
of the complex, nothing above the second dimension is 
required. 

The complexity of computing the link of w is directly 
related to the number of simplices containing w as a 
vertex. In fact, given a list of all simplices in TZ which 
include vertex w, print the simplex obtained by remov- 
ing vertex w from the simplex. This is the link, so given 
a listing of the simplices with vertex w, the link of w 
can be computed in time 0{k) where k is the number 
of adjacent simplices. Since k is obviously less than ro, 
the time to initially compute all the links is 0(mn). 

The time to compute Hi(Lk(w)) is 0(k 3 ), again with 
k equal to the number of simplices in the link of w. 
Since, k < m, the time taken for this initial step is 
0(nm 3 ). Thus, the entirety of the initialization step 
takes time 0(m 3 + ran + rim 3 ) = 0(nm 3 ). 

Failed Vertex 

In the worst case, every failed vertex is flagged and so 
the dS-G criterion must be recomputed each time. As 



seen in section 6.4 updating the matrix R and check- 
ing the dS-G criterion takes time 0(t 2 ) where t is the 
number of two simplices. 

What is interesting about the link condition is that it 
is easy to maintain the links of all the interior vertices. 
Let Lk((7, X) be the link of a in the simplicial complex 
X and let X w be the largest subcomplex of X without 
the vertex w. 

Lemma 4. For any vertices v,w in a simplicial complex 
X, 

Lk(v,X w ) = Lk(v,X)nX w . (2) 

Proof. If er € Lk(v,X w ) then obviously a E X w . Also, 
this implies that the simplex t =< a,v >G X w . Since 
a < t and v ^ a, we must have that a £ Lk(w,A), so 
Lk(u, X w ) C Lk(u, X) n X w . 

Let a <E Lk(w, X) n X w . Since it is in Lk(w, X), the 
simplex, < a, v >€ X. As a is also in X w , w £ a, so 



Combining this with the initializing step, the whole al- 
gorithm takes at worst time 0(m 3 n 2 ), so is polynomial 
in the number of simplices. 

8 Conclusions and Possible Ex- 
tensions 

In this paper, we have extended the problem posed by 
de Silva and Ghrist in [9] by assuming that sensors have 
a probability of failure, and asking for the probability 
of failure of the dS-G criterion for coverage. We de- 
termined that the generalized version of the problem is 
#P-complete, and thus it is unlikely that there is an 
algorithm to answer this question in general which runs 
in a reasonable amount of time. Finally, we provided a 
deterministic algorithm which does work in the case of 
a small set of sensors, and a method to predict failure 
when the system is larger but is being monitored. 

The obvious immediate extension of our work is 
to determine whether the the version of the problem 
posed by de Silva and Ghrist in 
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which allows for 
higher dimensions and looser boundary conditions, is 
also amenable to an application of probability of fail- 
ure. We conjecture that this extended problem will also 
be NP-hard. 

In the long term, we would like to see more applica- 
tions of computational topology to the design and anal- 
ysis of sensor networks. Since we can make such strong 
conclusions with such weak assumptions on the capa- 
bilities of the sensors, we expect that such applications 
are abundant. 

References 

[1] PK Agarwal, H Edelsbrunner, J Harer, and 
Y Wang. Extreme elevation on a 2-manifold. 
Discrete & Computational Geometry, 36:553-572, 
2006. 

[2] YA Ban, H Edelsbrunner, and J Rudolph. Inter- 
face surfaces for protein- protein complexes. In Pro- 
ceedings of the eighth annual international confer- 
ence on Resaerch in computational molecular biol- 



11 



ogy, RECOMB '04, pages 205-212, New York, NY, 
USA, 2004. ACM. 

[3] G Carlsson. Topology and data. Bulletin of 
the American Mathematical Society, 46(2):255- 
308, January 2009. 

[4] G Carlsson, T Ishkhanov, V de Silva, and 
A Zomorodian. On the local behavior of spaces of 
natural images. International Journal of Computer 
Vision, 76:1-12, 2008. 10.1007/sll263-007-0056-x. 

[5] EW Chambers, V dc Silva, J Erickson, and 
R Christ. Vietoris-Rips Complexes of Planar 
Point Sets. Discrete & Computational Geometry, 
44(l):75-90, July 2009. 

[6] F Chazal, LJ Guibas, SY Oudot, and P Skraba. 
Analysis of scalar fields over point cloud data. In 
Proceedings of the twentieth Annual ACM-SIAM 
Symposium on Discrete Algorithms, SODA '09, 
pages 1021-1030, Philadelphia, PA, USA, 2009. So- 
ciety for Industrial and Applied Mathematics. 

[7] D Cohen-Steiner, H Edelsbrunner, and D Morozov. 
Vines and vineyards by updating persistence in lin- 
ear time. Proceedings of the twenty-second annual 
symposium on Computational geometry - SCG '06, 
page 119, 2006. 

[8] CJ Colbourn. The Combinatorics of Network Re- 
liability. Oxford University Press, 1987. 

[9] V de Silva and R Christ. Coordinate- free Coverage 
in Sensor Networks with Controlled Boundaries via 
Homology. The International Journal of Robotics 
Research, 25(12):1205-1222, December 2006. 

[10] V de Silva and R Christ. Coverage in sensor net- 
works via persistent homology. Algebraic & Geo- 
metric Topology, 7:339-358, April 2007. 

[11] H Edelsbrunner and J Harer. The persistent morse 
complex segmentation of a 3-manifold. In Nadia 
Magnenat-Thalmann, editor, Modelling the Phys- 
iological Human, volume 5903 of Lecture Notes in 
Computer Science, pages 36-50. Springer Berlin / 
Heidelberg, 2009. 

[12] H Edelsbrunner and J Harer. Computational Topol- 
ogy: An Introduction. American Mathematical So- 
ciety, 2010. 

[13] MR Garey and DS Johnson. Computers and 
Intractibility: A Guide to the Theory of NP- 
Completeness. W.H. Freeman and Co., 1979. 

[14] R Ghrist and A Muhammad. Coverage and 
hole-detection in sensor networks via homology. 
IPSN 2005. Fourth International Symposium on 
Information Processing in Sensor Networks, 2005., 
(1):254 260, 2005. 



[15] JJ Headd, YEA Ban, P Brown, H Edelsbrunner, 
M Vaidya, and J Rudolph. Protein-protein in- 
terfaces: properties, preferences, and projections. 
Journal of Proteome Research, 6(7): 2576-2586, 
2007. PMID: 17542628. 

[16] A Tahbaz-Salehi and A Jadbabaic. Distributed 
coverage verification in sensor networks without 
location information. Automatic Control, IEEE 
Transactions on, 55(8):1837 -1849, aug. 2010. 

[17] LG Valiant. The complexity of computing the per- 
manent. Theoretical Computer Science, 8(2): 189 - 
201, 1979. 

[18] J Yick, B Mukherjee, and D Ghosal. Wire- 
less sensor network survey. Computer Networks, 
52(12):2292 - 2330, 2008. 



12 



